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I. INTRODUCTION 



In recent years, the study of quantum mechanical systems with a position-dependent 
mass (PDM) have attracted a lot of interest due to their relevance in describing the physics 
of many microstructures of current interests, such as compositionally graded crystals 
semiconductor heterostructure 0, quantum dots jjjj], 3 He clusters (^^metal clusters jBf etc. 
The concept of PDM comes from the effective-mass approximation p, @] which is a useful 
tool for studying the motion of carrier electrons in pure crystals and also for the virtual- 
crystal approximation in the treatment of homogeneous alloys (where the actual potential is 
approximated by a periodic potential) as well as in graded mixed semiconductors (where the 
potential is not periodic). Recent interest in this field stems from extraordinary development 
in crystal-growth techniques like molecular beam epitaxy, which allow the production of 
nonuniform semiconductor specimen with abrupt hetero junctions In these mesoscopic 
materials, the effective mass of the charge carrier are position dependent. Consequently, the 
study of the position-dependent mass Schrodinger equation (PDMSE) becomes relevant for 
deeper understanding of the non-trivial quantum effects observed on these nanostructures. It 
has also been found that such equations appear in many different areas. For example, it has 
been shown that constant mass Schrodinger equations in curved space and those based on 
deformed commutation relations can be interpreted in terms of PDMSE [9]. The PDM also 



appear in nonlinear oscillator [10J, [11] and PT-symmetric cubic anharmonic oscillator [12] . 



The most general form of the PDM Hamiltonian proposed by von Roos [13] is defined by 

H — — i (^m(q) a ^-m(qf^-m(qy + m(qf ^m{qf ^m(q)°^j + V{q), (1.1) 

where the ambiguity parameters a, (3, 7 are related by a + /? + 7 = —1. The above 
Hamiltonian always has the following form: 

H = ~ 1 2^^ + 1 2^^Tq + U{ " q)l (1 ' 2) 
where the effective potential U(q) is given by 

„(, ) = K(,)-(a + 7 )^ + (or, + a + 7 )^. (1 . 3 ) 

It is quite natural that physical interests just described above have also enhanced the 
studies on exact solutions to PDMS E |14| -l30l] by employing various methods e.g . supersym- 



metric (SUSY) quantum mechanics [3JJ and point canonical transformation [32| to mention 
a few. Later, PDM quantum systems were successfully formulated in the framework of 



A/"-fold SUSY in Ref. 33], which has provided until now the most general tool for con- 
structing a PDM system which admits exact solutions because of its equivalence to weak 
quasi-solvability. To avoid confusion, we here note that A/"-fold SUSY is different from non- 
linear SUSY which has been long employed since the work by Samuel and Wess 34| in 
1983 to indicate the nonlinearly realized SUSY originated from the work by Akulov and 
Volkov [35] in 1972. For a review of TV-fold SUSY see Ref. [36], while for recent works on 
nonlinear SUSY see, e.g., Ref. [37J and references cited therein. 

Very recently, new classes of exactly solvable PDM quantum systems whose eigenf unctions 



are expressible in terms of so-called Xi polynomials were constructed in Ref. [28|. The new 
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findings of X n polynomials (n > 1) were associated with the more fundamental mathematical 



concept of exceptional polynomial subspaces of codimension n introduced in Refs. [381-140 



whose origin can be traced back to the pioneering work on the classification of monomial 



spaces preserved by second-order linear differential operators [41 . 

The purpose of the present paper is two-fold. The first one is to bring the purely mathe- 
matical concept of exceptional polynomial subspaces into more physical settings by allowing 
the position dependence of mass (in a spirit similar to Ref. j28[) in the framework of M- fold 
SUSY. In the constant-mass case, form of potentials related to exceptional polynomial sys- 
tems is very limited. Thus, we can enlarge the physical applicability of the mathematical 
concept by introducing PDM to quantum systems. On the other hand, the framework of 
A/'-fold SUSY enables us to talk about the physical phenomenon of dynamical A/'-fold SUSY 
breaking. The second purpose is actually to examine effect of PDM profiles on dynamical 
breaking of A/'-fold SUSY. In this respect, it is rather surprising that there have been few 
papers, like Ref. jllj], where broken as well as unbroken SUSY is described in PDM back- 
grounds depending on the mass profiles. One of the main reasons would be that SUSY has 
been mostly used just as a technique to obtain exact solutions. The true significance of 
the Witten's SUSY quantum mechanics j3l[, however, rather resides in the nonperturbative 
aspects of dynamical SUSY breaking. Hence, one of our main purposes is, in other words, 
to examine change of nonperturbative nature of quantum systems caused by variations of 
mass profiles in view of dynamical AAfold SUSY breaking. 

The paper is organized as follows. In Section [TT1 we provide a self-contained review of 
A/'-fold SUSY in a PDM background, especially for those who are not familiar with the 
subject. We also summarize mathematical structure of type B and type X 2 A/'-fold SUSY. 
In Section IIII[ we construct several A/'-fold SUSY PDM quantum systems and examine 
dynamical A/'-fold SUSY breaking in different PDM backgrounds. The first three models of 
type B A/'-fold SUSY have rational, trigonometric, and exponential potentials in the constant 
mass case. We show in particular that the models whose bound state eigenfunctions were 
shown to be expressed in terms of X x polynomials in Ref. [42j for the constant mass case and 
in Ref. (281 for the PDM cases can be obtained as type B systems. The last three models of 
type X2 A/'-fold SUSY have rational, hyperbolic, and exponential potentials in the constant 
mass case. For both types of A/'-fold SUSY, we find that the rational potentials have steady 
A/'-fold SUSY against variation of mass profile while all the other types of potentials can 
receive effect of PDM on their dynamical breaking of A/'-fold SUSY. Finally, we summarize 
the results and discuss their implications and prospects in Section [TV] 



II. REVIEW OF TV-FOLD SUPERSYMMETRY IN A PDM BACKGROUND 

An A/^-fold SUSY one-body quantum mechanical system with PDM is composed of a pair 
of PDM Hamiltonians 

R ± = 1 d 2 m'{q) d . 

2m{q) dq 2 2m{q) 2 dq W ' K ' J 



and an A/"th-order linear differential operator 



if ^ m . , d 



^-W^^+E^ML, (2-2) 
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which satisfy the following intertwining relations 



P^H- = H+P N , P+H+ = H-P+. (2.3) 



In the above, P£ is the transposition [43J of Pj^ given by 

Pi = (Pjvf = (-^fm(q)-^ + £ (-^^r 1 ^. (2.4) 

Actually, the two relations in (I2.3P are not independent; the first implies the second and vice 
versa, since the PDM Hamiltonians (12. ip are invariant under the transposition (H ± ) T = H ± . 
One of the significant consequences of the intertwining relations ( 12.31) is weak quasi- 



sovability, that is, H ± preserves a finite-dimensional linear space spanned by the kernel 
of the operator 

ff^cV^, V±=kerP±. (2.5) 
Each space Vjw- is called a solvable sector of H ± . Except for the M = 2 case (cf., Refs. [3 



441]). virtually all the A^-fold SUSY systems so far found admit analytic expression of Vjy in 



closed form, and thus are quasi-solvable. In addition, it sometimes happens when either H 
or H + does not depend essentially on Af and preserves an infinite flag of the solvable sectors 

Vr /+ C V 2 /+ C-CV/c-. (2.6) 

In this case, it is said to be solvable, which is a necessary condition for exact solvability. We 
note that H~ and H + are usually simultaneously solvable due to the intertwining relations 
(12~3|). 

A set of an A/"-fold SUSY system H ± and P^ provides a representation of A^-fold super- 
algebra defined by 

[Q%, H] = {Q± , Q±} = 0, {Q^, Q + M ) = 2*P M (H), (2.7) 

where Pn(x) is a monic polynomial of degree M in x. Indeed, it is realized by defining H 
and Q^f as 

H = H~r^ + + H + ^ + r, Qtf = P^ + , Qjr = Pfrl>-, (2.8) 

where tp is a pair of fermionic variables satisfying {?/' ± ,'^ ± } = and = 1. It is 

easy to check that the above H and satisfy the first part of algebra (12. 7p . In particular, 
the intertwining relations in (12. 3p guarantee the commutativity of H and Q%- Regarding 
the second part of algebra, the monic polynomial Pj\f is given, in the above representation, 
by HI, S3 



P M (H) = det^H-H ± \ v ^j, (2.9) 

namely, the characteristic polynomial for H ± restricted to the solvable sectors V^. 

Whether A/"-fold SUSY of the system under consideration is dynamically broken is deter- 
mined by a property of the solvable sectors Vjy since they characterize AA-fold SUSY states, 
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namely, states annihilated by the pair of TV-fold supercharges Q^. Let |0) and |1) be the 
fermionic vacuum and the one fermion state, respectively, which satisfy 

Vlo) = o, |i) = ^ + |o). 



(2.10) 



Then, superstates |\l/ ) = \& (g)|0) and \^q) = ^o(?)|l)> respectively, are annihilated by 
both of Q% 



Q±\%) = 0, Q±-|* + >=0 



(2-11) 



if and only if (?) £ V^- and ^o(q) G V^, respectively. However, such states do not 
necessarily satisfy physical requirements. Suppose S C C is a domain where both of the 
Hamiltonians H^ 1 have no singularities and are thus naturally defined, and 5(5') is a linear 
space of complex functions in which both of H^ 1 act. In a usual physical application, the 
domain 5* is the real line Mora real half-line R + = (0, oo), and the linear space 5 is a Hilbert 



space L 2 , so that 5(5) = £ 2 (R), or L 2 



In the latter cases, the physical requirement is 



the normalizability (square integrability) on 5. Then, there exists physical (normalizable) 
TV-fold SUSY states \%) and/or which satisfies fl2TTTj) if V^(S) C L 2 (S) and/or 

V^(5) C L 2 (S), in other words, if H~ and/or H + is quasi-exactly solvable. If there are no 
such physical TV- fold SUSY states in the Hilbert space L 2 (5) exists then TV- fold SUSY of 
the system is said to be dynamically broken. It was first shown correctly in Ref. |45(] that the 
generalized Witten index characterizes TV-fold SUSY breaking, which corrected the wrong 
statement made earlier in Ref. [46| . 

For TV > 1, we can have an intriguing situation where not the whole of, but a subspace 
of the solvable sectors V^-(S) and/or V^(5) belong to the Hilbert space L 2 (5). In this case, 
TV-fold SUSY of the system is said to be partially broken. Partial breaking of TV-fold SUSY 
was first discovered in Ref. (47|. We note that it is different in nature from the partial 
breaking of (nonlinear) SUSY 0,111. 

Construction of an TV- fold SUSY system is in general quite difficult, especially for a larger 
value of TV, since the intertwining relations (12.31) compose of coupled nonlinear differential 
equations for U ± {q) and w^(q) {k = 0, . . . ,TV — 1). For the direct calculations of inter- 
twining relations in a PDM background in the cases of TV = 1 and 2, see Ref. 24j. To 
circumvent the difficulty, a systematic algorithm for constructing an TV-fold SUSY system 
was developed in Ref. 47 [for constant-mass quantum mechanics and was later generalized 
to PDM systems in Ref. 33j. The significant feature which is common in both constant-mass 
and PDM systems is that an TV- dimensional linear space of functions 



V 



(2.12) 



preserved by a second-order linear differential operator H~ can determine whole of an TV-fold 
SUSY system. Indeed, we can construct a pair of TVth-order linear differential operators Pj^ 

and another TV- dimensional vector space Vjjr such that Vj^ = ker P^. 
that a pair of second-order linear differential operators given by 



Then, we can show 



H 



± 



d 2 

az 2 



TV - 2 



A'(z)±Q(z) 



C(z) 



:i±d 



TV-l 



Q'(z) 



-A\z)w%Uz) - A{z)w^{z) 



,7,W' 



(2.13) 
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is weakly quasi-solvable with respect to the spaces VV, namely, H^V^ c ^Af- 



by 



With the choice of the change of variable z = z(q) and the gauge potential Wj^ determined 



z'(q) 2 = 2m(q)A(z) 



In 



we can obtain an A/"-fold SUSY system by 



e -^ H ± e w^ 



z=z(q) 



P ± 



z=z{q) 



(2.15) 



With the change of variable and the gauge transformation, both of H ± get the form of PDM 
Hamiltonian (12. ip and their effective potentials U ± (q) are given by 



1 



2m(q) 

-a±i 



dW 



dq 



d 2 W 



dg ; 



v | m , (g)dny 
m(g) dg 



<7(s(g)) 

„7,M' 



(2.16) 



*(<?) 



The solvable sectors Vt? of i? 1 * 1 are evidently given by 



(2.17) 



In principle, we can construct a pair of A/"-fold SUSY PDM Hamiltonians if 1 * 1 and its solvable 
sectors Vjy by using the formulas ( I2.16P and (I2.17p . However, there is an easier way to obtain 
such a system when we have already had an ordinary A/"-fold SUSY constant-mass quantum 
system at hand. Suppose the latter system is such that its pair of potentials V^ ± (q), its 
gauge potentials VV^' )± (g), its solvable sectors V$ [q] are all known. Then, an A/"-fold 
SUSY PDM system having a pair of effective potentials U ± {q) ) gauge potentials VV^(g), 
and solvable sectors V^[g] can be constructed immediately via the following prescription: 



u±{ ) = V (o)± (u(q)) + m "^) _ 7m '(l) 2 
W V [U[q)) + 8m(qy 32m(g)3- 



V^( ff ) = --ln|m(g)|+W v 

V^[q\=m{q) 1 "V l ff ± [u{q)i 
where the function u(q) is given by 



(0)±, 



«(?))> 



(q) = / dqy/m(q). 



(2.18a) 
(2.18b) 
(2.18c) 

(2.19) 



Actually, the above relations are consistent with the formulas obtained by the point canonical 
transformation, see, e.g., equations (2.7 ) and (2.8) in Ref. [22|, equation (7) of 21] and 
equations (10), (13), and (14) in Ref. [19|. The above relations (I2.18P have also been verified 
in Ref. [33] where type A A/"-fold SUSY has been constructed in PDM background. One 
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of the most salient features unveiled by the algorithmic construction is that both constant- 
mass and PDM quantum systems with A/"-fold SUSY have totally the same structure in 
the gauged z-space. That is, the functional forms of the gauged operators such as and 
H^ 1 given by (12.131) are identical in both the cases. It means in particular that the starting 
vector space Vjj- determines all in the algorithm regardless of whether mass is constant or 
not. Hence, different types of A/"-fold SUSY are characterized by different types of vector 
spaces Vy and vice versa. Until now, four different types have been discovered, namely, type 
A H Ej, type B [51], type C 0, and type X 2 0. We note that almost all the models 



having essentially the same symmetry as A/"-fold SUSY but called with other terminologies 
in the literature, such as Poschl-Teller and Lame potentials, are actually particular cases of 
type A A/"-fold SUSY. In this article, we focus on constructing PDM quantum systems with 
type B and type X 2 A/"-fold SUSY since the other types (type A and type C) are not related 
to exceptional polynomial subspaces. In what follows, we shall review the general structure 
of these two types of A/"-fold SUSY. 



A. Type B AA-fold Super symmetry 



Type B A/"-fold SUSY was first discovered in Ref. [5JJ and was found to be associated 
with the following monomial space 

V^ = V^:=<W..,^- 2 ,^}, (2-20) 



called type B, which was considered in Ref. [41] in the context of the classification of mono- 
mial spaces preserved by second-order linear ordinary differential operators. Applying the 
algorithm to the type B monomial space, we obtain j36[ the gauged A^-fold supercharge 
components 

and the functions which characterize the gauged Hamiltonians (12 . 1 3[) are given by 

A(z) = a^z 4 + a^z 3 + a 2 z 2 + a\Z + ao, (2.22) 
2Q(z) = -Ma z z 2 + 2b lZ - Na x , (2.23) 
C(z) = Af(Af - 3)a 4 z 2 + A/VV - 2)a 3 z + c , (2.24) 

and w <L ^\ 1 (z) = —z^ 1 . The other linear space preserved by H + is given by 

V+- = z- l (l,z 2 ,...,z"). (2.25) 

We note that both the monomial spac es (I2.20p and (I2.25P are actually exceptional polynomial 
subspaces of codimension 1, see Ref. 40j. We can easily check that the type B Hamiltonian 
H + preserves an infinite flag of the following spaces 

y+ e -< c V 2 + e~< C ■ ■ ■ C V>~< C ■ ■ • , (2.26) 

where Vjjr and VV^ are given by (12.251) and (12.141) . respectively, and thus H + is solvable if 
and only if = = 0. On the other hand, the partner type B Hamiltonian H~ does 
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not appear to be solvable for any parameter value since the type B monomial space ( I2.20p 

does not constitutes an infinite flag due to the fact that vffi (f. Vjy!L for all M = 1, 2, 

However, it turns out [36] that, when a 3 = a 4 = and H + gets solvable, the partner 
Hamiltonian H~ does preserve an infinite flag of linear spaces given by 



where VV^ is g: 



iven by (12.141) and "Pj^ is the type A monomial space defined by 
V^ = (l,z,...,z^). 



[2.27) 



(2.28) 



That is, H and H + can be solvable simultaneously. In this paper, all the type B models 
we will consider later satisfy the solvability condition = 0. Thus, all the pairs of 

type B Hamiltonians if ± preserve the infinite-dimensional solvable sectors V given by 



V- = <l,z(g),z(g) 2 ,...)e- w ^) 
V + = (l,z(q) 2 ,z(q)\...)z(q) 



(2.29) 



An interesting consequence of the fact that H and H + preserve different types of infinite 
flag of spaces in the solvable case is that the eigenfunctions of H~ are expressed in terms of 
a classical polynomial system while those of H + are in terms of an X\ polynomial system. 
It is exactly the underlying reason why some of the Hamiltonians whose eigenfunctions 
are expressed in terms of the X\ Laguerre or Jacobi polynomials were obtained by those 
whose bound state eigenfunctions are expressed in terms of the classical Laguerre or Jacobi 
polynomials using an intertwining or SUSY techniques in Refs. [2?], 53, 54]. 



B. Type X2 AA-fold Super symmetry 



Type X2 A/"-fold SUSY constructed in Ref. 52j is associated with the following exceptional 
polynomial subspace of codimension 2 

= {<Pi(z\ <*),..., 0tf(z;a)), (2.30) 

where y5 n (z; a) is a polynomial of degree n + 1 in z with a parameter 0, 1) defined by 

<p n {z; a) = (a + n- 2)z n+1 + 2(a + n - l)(a - l)z n + (a + n)(a - l)az n-1 . (2.31) 



Applying the algorithm to the X 2 space (I2.30p . we obtain [5^] the gauged M- fold supercharge 
components 



f{z;a) 



AT-l 

n 



f(z; a + k + l) ( d f(z; a + k + 1] 



P+ 



z (q) 



f(z;a + Af) ^ f(z;a + k) \dz f(z;a + k + l) 
Y^(± + f'(z;a + Af-k)\ f(z;a+M-k) 



k=0 



dz f(z;a + J\f-k)Jf(z;a+J\f-k-l 



f(z;a) 



(2.32) 



f(z;a+J\f)' 
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where n^o 1 '■= A/V-i ■ ■ ■ AiA , and the functions f(z; a) and w^f^^z) are given by 

f(z;a) = z 2 + 2{a-l)z + (a-l)a, (2.33) 

v>i,Uz) = -CAT " - f{z . a + U y ( 2 - 34 ) 

The most general forms of the functions A(z), Q(z), and C(z) appeared in depend on 

four parameters Oj (i = 1, . . . , 4), but in this paper we only consider models with = a 3 = 0. 
In the latter case, they read as 

A(z) = a 2 z 2 + aiz + (a - l)(a + N - l)a 2 , (2.35) 

Q(z) = — a 2 z 2 — (3a 2 + a i)z — (a — l)(3a + 3Af — 7)a 2 
+ 4(a-l)0 W 
2 /(z;a) 
4(a-l)D(z) 

C z = a 2 z + c - 1 ,, V » ( 2 - 37 ) 

J{z;a) 

where D(z) is given by 

D(z) = -[(2a +N- 3)a 2 - a t ]z - (a - l)(2a + jV - l)a 2 + aai. (2.38) 

For their most general forms, please refer to Ref. [52] . The other linear space preserved 
by H + is given by 

V£ = (xi(z; a + AT), ■ ■ • , X^; « + AO)/ 0; a + A/") -1 , (2.39) 

where Xn(z; a ) is a polynomial of degree n + 1 in z defined by 

X n (z] a) = (a — n)(a — n + l)z n+1 + 2(a - n — l)(a — n + l)(a - l)z n 

+ (a-n- l)(a-n)(a- l)a2 n_1 . (2.40) 

The solvable sectors Vjy of the constant-mass Hamiltonians H ± are 

V M = (<pi(z(q); a),..., m{z(q); a)) e~ w ^\ (2.41) 
v+ _ (xMq); a + AQ, . . . , x^(*(g); a + A/")) w + (g) 

" /(%); + AO ' ^ 

Finally, the type X 2 Hamiltonians H ± preserve the infinite flag of the spaces Vjy (A/" = 
1,2,.. .) and are simultaneously solvable if and only if a 2 = (a 3 = a 4 =)0. 



III. TYPE B AND TYPE X 2 TV-FOLD SUPERSYMMETRY FOR POSITION- 
DEPENDENT MASS 

In this section, we shall consider some models which belong to type B and type X 2 M- 
fold supersymmetry. In order to study effect of PDM in these models, we need to consider 
simultaneously the corresponding constant-mass type B and type X 2 models as well. In 
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particular, we shall address ourselves to the following question: Does position dependent 
mass have any effect on dynamical breaking of type B and type X 2 TV-fold SUSY? By 
comparing the solvable sectors of both constant and position-dependent mass cases, we 
shall see below that the answer is in the affirmative in some cases for particular choices 
of physically interesting mass functions. In order to explore in detail the impact of mass 
functions on symmetry breaking or restoration, it will be appropriate to consider more 
than one mass function in a few examples. Also it will be shown that the bound state 
wavefunctions of one of the partner potentials obtained in type B TV-fold SUSY are associated 
with exceptional Xi Laguerre and Jacobi polynomials while those of the other partner are 
associated with classical Laguerre and Jacobi polynomials. 



A. Effects of PDM on Dynamical Symmetry Breaking of Type B TV-fold SUSY 

Here we shall consider three examples of type B TV-fold SUSY corresponding to 
three different choices of A{z). In each of the examples, we first show the results in 
the constant mass case, followed by the corresponding results in the PDM case. As we 
referred to before, all the type B models constructed below satisfy the solvability condi- 
tion and thus their solvable sectors in the constant-mass case are given by (12.291) . 

Example 3.1. A(z) = k(z - z ) (k ^ 0) 
Potentials: 

^(,)4^ + 4( ^-^ )2 - fc % ^ + ^ (3.1) 



V ^+(a) = ^a 2 + 4 W ~ k + 2k - 8kZo + V 
yq) 8 q 8k 2 g 2 kg 2 + 2z {kg 2 + 2z ) 2 °' 



(3.2) 



where Vo is an irrelevant constant given by 

~2k A' 

Solvable sectors: 



Vo = r T7 — it- 



V (0) " = (1, z(q), z(g)\ . . .) q ^ob 1 -2Mk + k)/(2k) e b iq y^ (3>3) 
V (0)+ = (1, z(g) 2 , z(g) 3 , . . >(g)~ V (2 ^ fc)/(2fc) e- 6ig2/4 . (3.4) 

We assume k > and zq > so that the pair of potentials V ± (g) has no singularities except 
for at q = 0. Thus, the system is naturally defined in L 2 (1R + ), IR + = (0,oo). In the latter 
Hilbert space, V (0) ~(M+) C L 2 (R + ) if and only if 

61 < and zoh > (TV - l)k, (3.5) 

which cannot be satisfied by any bi E R. On the other hand, V (0)+ (M + ) C L 2 (R + ) if and 
only if 

bi > and k > z bi. (3.6) 
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Hence, TV-fold SUSY of the system is unbroken if and only if < b\ < k/zo on the constant- 
mass background. 

Now, the relevant expressions for partner potentials, gauge potentials and corresponding 
solvable sectors of type B PDM systems can be obtained using Eqs. (I3.ip - fl3.4p and relations 
f)2.18p . Since our main objective in this section is to study effect of mass function on 
dynamical breaking of A^-fold SUSY, we give below only the solvable sectors for an 
arbitrary mass function m(q): 

V- = {l,z{u{q)),z{u{q))\..)m{q) l l\{qf 2z «^ (3.7) 
V+ = (1, z(u(q)) 2 , z(u(q)) 3 , . . ^(^(g))-^^) 1 / 4 ^^)-^" 61 ^/^^-' 1 ^ 2 / 4 , (3.8) 

where u(q) is given by fl2.19p . At this point, we are in a position to choose a particular mass 
function. Let the mass function be 

m{q)=e- bq , b>0, gG (-00,00), (3.9) 

which was considered in Ref. j28[ where the PDM potentials were associated with X\- 
Laguerre polynomials. This exponentially behaved mass function has been often used in the 
study of confined energy states for carriers in semiconductor quantum well [l9|, 28]. It has 
also been used to compute transmission probabilities for scattering in abrupt heterostruc- 
tures [25| which may be useful in the design of semiconductor devices j55j. For the mass 
function, the change of variable is given by 

u(q) = ~\z- hql \ (3.10) 
and the pair of potentials U ± (q) reads from (I2.18a|) as 

<T <*) = + -*V + + Vo , (3,1) 

U {q) ~2V e + 8P 6 + 2ke- b « + z b* (2fce-^ + z & 2 ) 2 °' 1 ' 

respectively. It is worth mentioning here that the potential U + (q) given in fl3.1ip is iden- 
tical with the potential V e ff(q) associated with exceptional Xi Laguerre polynomials [e.g., 
Eq. (12) of Ref. (Hj]], if one takes k = 1/2, 61 = b 2 /2, and z = a/b 2 . On the other hand, 
for the same choices of parameters the other potential U~ (q) coincides with the potential 
[after making a translation a — > a — Af] previously obtained in Ref. [27| corresponding to 
classical Laguerre polynomials. 

The solvable sectors of the potentials (13. lip and (I3.12p are respectively given by 



V- = <l,e-^ + ^ ,(e- fe9 + ^o) 2 ,...> 



x exp 



(3.13) 



<l,(e-^ + ^ ) 2 ,(e- b9 + ^o) 3 ,...> 



x e 



-bq 



+ zq) exp 



k 2 q & 2 



(3.14) 
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where z = z b 2 /(2k). Here the potentials have no singularities in the finite part of the real 
line, so the domain is R. Since b > 0, so V~(R) C L 2 (R) if and only if b\ < 0. On the 
other hand, V+(R) C L 2 (R) if and only if b x > 0. Hence, the Af-fold SUSY of the PDM 
system is unbroken unless b\ — 0. Comparing the solvable sectors of both the constant and 
position-dependent mass scenarios, it can be observed that it is not possible to break TV-fold 
SUSY dynamically for the particular choice of mass function m(q) = e~ bq . In addition, 
we have checked that many physically interesting mass functions also have no effect on 
symmetry breaking. 

Example 3.2. A(z) = a 2 [l - (z - z ) 2 ]/2 (a > 0) 
Potentials: 

T ,roW x W -Af 2 a 4 )z smaq (2b x - Ma 2 ) 2 z 2 + (2b x + Afa 2 ) 2 - a 4 2 

V { °>~ (q) = A 9 — ^~ + — - \ tan 2 aq 

Aa z cos z aq 8cr 

+ ^ + (3.15) 

(0)+ (2&J - Afa 2 ) 2 z smaq {2b x - Afa 2 ) 2 {z 2 + 1) - a\ 2 

v (?) = 5 1 ^ tan aa 

Aa z cos z aq 8cr 

a2z o « 2 (^o 2 ~ jj b\J\f 

+ : t- 1 \2 T~ + V °> ( 3 - 16 ^ 

sin aq + z (sin aq + Zq) 2 



where Vq is an irrelevant constant given by 

h a 2 (Af 2 - 7 
Af + 12 8a 2 



_ . h a 2 (N" 2 - 7) (2Mo - Afz a 2 ) 2 

Vq = -r? H -— 1 -rt. 



Solvable sectors: 



i I .y-i / 1 + sin ag 



(26 1 -ATa / )z 



4 a 



V(°)- = <1, z(g), z(g) 2 , . . .)| cosag|^+^ , (3.17) 

v ' \ 1 — sm aq J 

(2b 1 -ATa 2 )z 

V<°> + = (l,z( q ) 2 ,z( q f,...)z( q r\co S a q r^^ (1±^L) = . (3.18) 
x ' \ 1 — sm aq J 

It is worth mentioning here that the potential V^°' + ( q ) coincides with the potential whose 
bound state wave functions are given in terms of exceptional X\ Jacobi polynomial |42| for 
a = 1, &! = B +N/2 L zo = -(2 A - 1)/(2B) whereas potential U (0)_ (g) coincides with 
the Scarf I potential [54| [after making an change B — )■ B + A/"] whose bound state wave 
functions are given in terms of classical Jacobi polynomials. 

We choose here a domain of the system as S = and assume z > 1 so that the 

pair of potentials V^ ± (q) has no singularities except for at the boundary OS = { — ^}. 
Thus, the Hilbert space for the system is L 2 (S). Then, V ( -°- ) ~(S') C L 2 (S) if and only if 

&! AT - 1 (2&i - A/a 2 )z 1 
^ + 2 2^ > ~2' 



that is, 



Afa 2 z - 1 , Afa 2 z + 1 , n , Q . n . 

^^TT <6l< ^^T for -o>l- (3.19) 
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Similarly, V (0)+ (S) C L 2 (S) if and only if 



b x Af-1 (2b 1 -Afa 2 )z 1 



that is, 



\ r 2 

h > and z > 1. (3.20) 

Hence, A/"-fold SUSY of the system is broken for the constant mass case if and only if Zq > 1 
and 

A/a 2 z - 1 A/V 2: + 1 

61 -^^n or i^r (3 - 21) 

In a PDM case, the solvable sectors V 1 * 1 of the type B PDM A^-fold SUSY partner Hamilto- 
nians H for an arbitrary mass function m(g) are deformed according to fl2.18c[) as 

V" = (1, z(u(q)), z(u(q)) 2 , . . .)m(q)* 

_(26 1 -AAa 2 )z 

/ m 6 i 4- K=i / 1 + sin auiq) \ 

x cos au(q) ^ + — : f± , 3.22 

\1 — smau{q) J 



V+ = (l,z( U (q)) 2 ,z(u(q))\...)m(q) l ± 



'l 1 AM 

coscra(g)| ^ 2 / 1 + sina-u(g) \ ^ 



x J : ^ ; p( . (3.23) 

sin au(q) + z \ 1 — sm att(g) / 

where w(g) is given by f)2.19p . In this case, the choice of mass function and the corresponding 
change of variable are given by 

2 2 

771(g) = -e~ 2q , u(q) = Erf q, q G (-00,00). (3.24) 

7T 

Consequently, the partner potentials U ± {q) read as 

_ {4b 2 - N 2 a A )z sin(aErf q) (3g 2 + l)^ 2 b x N 
^' ~ 4^2 cos 2 (a Erf g) 4 + ~Y~ + 

+ (^-^'W + (^+W)'-»' tm3(nEtf9) , (3.25) 



f/+(g) 



(26i - ATa 2 ) 2 z sin(a Erf g) (3g 2 + l)7re 2g2 b x N 



4a 2 



cos 2 (a Erf g) 



+ ^0 



+ 



a 2 z 



a 2 (z 2 



1) 



sin(a Erf g) + z [sin(a Erf g) + z ] 2 

(26 1 -A^a 2 ) 2 (z 2 + l)-a 4 2/ „ p N 
+ — " tan ( a Erf ?) ■ 



(3.26) 
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The solvable sectors of the potentials (I3.25P and (13.261) are given by 



V- = <l^Kg))^Kg)) 2 ,...>e-" 2 / 4 

H ! m-i f 1 + sin(a Erf q) 



x | cos(aErf q)\~> 



1 — sin(a Erf q) 



An'A 



(3.27) 



[l,z(u(q))\z(u(q)T,... 

bi i JV— 1 

| cos(aErf q)\ ~*~ 



sin(a Erf q) + zq \ 1 — sin(a Erf g) 



-9 2 /4 



1 + sin(a Erf g) 



(2b 1 -A"a^)z 
4a? 



(3.28) 



The potentials U ± (q) as well as the mass function are well behaved in q G (—00,00). So, 
we can take the domain as the whole real line R. Since Erf q — > ±1 as q — > ±00, so both 
the solvable sectors V ± (IR) belong to L 2 (IR), irrespective of the parameter values of bi and 
zq. Hence, it manifests unbroken SUSY. So, in this case position-dependent mass affects 
the symmetry breaking scenario. But the mass profile m(g) = sech 2 aq, q E (—00, 00) has 
no effect on dynamical breaking of A/"-fold SUSY which can be observed by considering the 
leading behavior of the solvable sectors (I3.22p and ( I3.23p . We have found that same is true 
for many other mass functions. 

Also associated to this mass profile, one of the partner potentials given in equation (I3.29p 
is identical with the V e ff(q) whose bound state wave functions are given by exceptional 
X\ Jacobi polynomials [e.g., Eq. (18) of Ref. |28|], for the choice of parameters b\ = (a — 
(3 + J\f)a 2 /2, zq = (a + (3) /(a — 0). The simplified form of the other partner potential 
U~{q) matches with the potential previously obtained in [2?| corresponding to classical 
Jacobi polynomials. It is worth mentioning that this mass profile m(q) = sech 2 aq has 
been previously used in PDM Hamiltonians of BenDaniel-Duke 5f| and Zhu-Kroemer 57 
type and interesting connection was shown [5B| between the discrete eigenvalues of such 
Hamiltonians and the stationary 1-soliton and 2-soliton solutions of the Korteweg-de Vries 
(KdV) equation. 

For the latter choice of the mass function, the change of variable is given by u(q) = 
tan _1 (sinhag)/a and corresponding pair of potentials U ± {q) read as 



[2ji (z + 1) - Afa 2 z =F (TV - 2)q 2 ] [26i (z + 1) - Afa 2 z =F (A/" + 2)q 2 ] 2aq 

32a 2 6 
[26i(z - 1) - Ma 2 z Q ± (JV - 2)a 2 }[2b l (z - 1) - J\fa 2 z Q ± (A/~ + 2)q 2 ] _ 2aq 

+ 32^ 6 



+ 



1 ± 1 a 2 
2 z + l 



1 - 



2(^o - 2) 



4(^o - 1) 



z - 1 + (zq + l)e 2a( i {zq-1 + (z + l)e 2a 9) 2 



=F — — + Vo- 



4 



(3.29) 



Example 3.3. A(z) = (z - -2 ) 2 /2 
Potentials: 



V^(q) 



(2b 1+ M) 2 z 2 (Ab 2 -Af 2 )z 

8 e -h 4 e -t- vo, 



v{Q)+{q) = (2b 1+ N) 2 z 2 c _ 2Q + (2b 1+ Af) 2 z c _ q _ 



z e~ 



1 + z e-i) 2 



+ V . 



(3.30a) 
(3.30b) 
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where Vq is an irrelevant constant given by 



6? &1 AP+11 

v °-^2 + Jr + ^^~ R - 



Solvable sectors: 

V(°)- = <M(g)^(g) 2 ,...)exp 

V^ + = (l,z(q) 2 ,z(qf,...)z(q)- 1 exp 



{2b x +M)z Q 



N - 1 - 2b x 



(2&i+A/>o 



-q 



N - 1 + 26i 



-9 



(3.31a) 
(3.31b) 



We assume > so that the pair of potentials V^ ± (q) has no singularities in (— oo, oo). 
As we will show in what follows, the Af-fold SUSY in this case can be partially broken. To 
see this, we first introduce a pair of fc-dimensional subspaces of the solvable sectors 

as 



V, 



(0)- 



(l,z{q),...,z{q) k - 1 )exp 



(2b 1 +Af)z 



V, 



(0)- 



(l,z{q) 2 ,...,z{q) k )z{q) 1 exp 



(26i + A>o 



M - 1 - 26i 



A/" - 1 + 2&i 



Then, for a fixed 6 N, we have 



v. 



(o)+. 



CL 2 
CL 2 



-A/" < 26i < jV + 1 - 2Jfe, 
2A; - A" - 1 < 26i < -A". 



(3.32) 
(3.33) 



(3.34) 
(3.35) 



From these conditions, it is easy to observe that (R) C L 2 (M) if and only if —N/2 < 
bi < (A" + 1 - 2k) /2 for a k G N satisfying k < N + 1/2, while there is no fc G N which 
satisfy the condition (EM and thus V (0)+ (M) £ L 2 (M) V&i G R. Hence, the jV-fold SUSY 
in the constant-mass background is partially broken if there is a positive integer k < M for 
which the parameter b\ satisfies 



N , N+l 

— < &i < 

2 2 



2A; 



and fully broken otherwise. 



The solvable sectors V^ 1 of the corresponding PDM Hamiltonians H are written as 

V- = (l,z(u( q )),z(u( q )) 2 ,...)m(q) 1 / 4 

[ (2j i + A>p _ u(g) _ A~ - 1 - 26i 



X6XP 2 2 

i^k?)) 2 ,^)) 3 ,...}^))-^) 1 / 4 



x exp 



(26 1 +AQgb c _ t , 



(3.36a) 



(3.36b) 



and the potentials U ± (q) can be obtained using Eqs. (12.18a[) . (13.30b ) and ( 13.30b ) . We 
have checked the normalizability of the solvable sectors (I3.36P with the following two mass 
functions. 
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(i) m(q) = (1 — q 2 )~ 1 , q G (—1, 1) for which the change of variable is u(q) = sin -1 q. This 
mass profile has been used in Refs. (Tol . 11 ] while considering the effective-mass quantum 
nonlinear oscillator. This mass function has effect on dynamical symmetry breaking because 
it manifests broken SUSY [i.e., neither V~ nor V + belongs to L 2 (— 1, 1)], which is clear from 
the following expressions of V~ and V + : 



V 



[l, z{u{q)), z{u{q))\ . . .)— — _ 



x exp 



(26i + Af)z 



i _ 



l-26i . _! ' 
sin g 



MKg)) 2 ,zKg)) 3 ,...; 
(2&!+tV>o 

x exp 



1 



1 



g2)l/4(eB 

tV-U 



r 9 + zb) 

26i . _j ^ 
sin q 



2 2 
(ii) m(g) = 2e~ 2qZ /n for which the solvable sectors (I3.36P reduce to 
V- = (l,z(u(q)),z(u(q)) 2 ,...) 



x exp 



' 4 



(26i +N)zq 



Erfg 



AT - 1 - 26i 



;i ) «(u(9)) 2 ) ... > z(tt(ff)rMu(g))- 1 



x exp 



"4 + 



(26i+j\02b 



e -Erfg _ 



TV - 1 + 26i 



Erf g 



Erf g 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



From the above solvable sectors, we observe that both V~(R) and V + (M) belong to L 2 (R), 
irrespective of the parameter value b x , which means unbroken TV- fold SUSY. Hence, the 
mass function m(g) = 2e~ 2q /it affects dynamical breaking of the TV-fold SUSY. 

Hence, comparing the normalizability conditions in both the constant and position de- 
pendent mass cases, we conclude that both the mass functions change the behaviours of 
symmetry breaking. 



B. Effects of PDM on Dynamical Symmetry Breaking of Type X 2 TV- fold SUSY 



In this section, we examine three different models of type X 2 TV-fold SUSY characterized 
by different choices of the two parameters a x and a 2 ; a x 7^ and a 2 = for the first model, 
a x — and a 2 7^ for the second, and a x a 2 7^ for the third. The first two choices lead to 



the rational- and hyperbolic-type potential pairs already shown in Ref. [52|, while the last 
choice to an exponential-type potential pair which is new and has not been investigated in 
the literature. 

Example 3.4. A(z) = 2z [a x = 2]. 
Potentials: 



V^ + {q) 



q 2 Aa 2 - 1 

— H h 4 

2 8q 2 

<f 4(q + TV) 2 - 
~2 + 8^2 



r q 2 -a + l 4(a-l)g 2 



jV;a) 2 
a-TV+l 



-TV+Uo, 

_ 4(a + TV- l)g 2 
/(? 2 ;« + TV) f{q^a + UY 



+ V 



(3.41a) 
(3.41b) 
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where Vq = -V — a + 3 — Cq is an irrelevant constant. 
Solvable sectors: 



W = (fdQ 2 ; a),..., w(q 2 ; «)> 



? a+l/2 e -g 2 /2 



< + = <Xi(<? 2 ; a + AT), ... , XM (q 2 - a + Af)) 



q -a-Af+l/2 e q 2 /2 



(3.42a) 
(3.42b) 



/(gV + AO • 

In this case, the solvability condition a 2 {= 03 = a 4 ) = for type X 2 is satisfied and thus the 
corresponding const ant- mass Hamiltonians H^- ^ are simultaneously solvable. 

For a > 1, a natural choice for the domain of these potentials is a real half-line S = M.+ . 
On this domain R+, it is evident from (137321) that V$~{R+) C L 2 (R + ) and Vj? )_ (R + ) <£. 
L 2 (R + ). Therefore, it manifests unbroken A/"-fold SUSY of the system in the constant-mass 
background. 

According to (I2.18cl) . the solvable sectors Vv of the corresponding PDM Hamiltonians 
H^ 1 for an arbitrary mass function m(q) read as 



V M = (vMq) 2 ; a),---, ipM(u(q) 2 ] a)) 



m(qy/*u(q) a+1 / 2 e- u ^ 2 / 2 



K = <XiM<?) 2 ; a + AT), ... , XM(u(q) 2 ; a + Af)) 



f(u{q) 2 ;a) 

m(g) 1 / 4 u ( g )-«-^+ 1 /2 e «(3) 2 /2 



/(t*(g) 2 ;a+jy0 



(3.43) 
(3.44) 



where w(g) is given by (12~T9|) and the PDM potentials U^q) can be obtained using 
Eqs. f)2.18ap and (13.4ip . In this case, we have not been able to find out any realistic mass 
function which could break the AA-fold SUSY. In other words, we can say that the A/"-fold 
SUSY in this case is steady against many variations of mass functions [e.g., m(q) = e~ g , 
sech 2 q] . 



- l)(a+AT-l) > 0]. 



Example 3.5. A{z) = (z 2 + C 2 )/2, [a 2 = 1/2, C 2 = (« 
Potentials: 

C 2 Af - 1 

V (0) ~ (q) = ^- cosh 2 q H — C sinh q + V 

8 4 

+ \A(Af-l)(smhq + 4a 2 + A(Af -2)a-AT 2 -2AT + A} 

8 cosh q 



2(a 



C sinh q - a-Af+3 , X 2C sinh q - Af + 1 
- — 2(a — V 



/(Csinhg;a) 



f((smhq;a) 2 



(3.45) 



C 2 SAT - 1 
V (0)+ (q) = ^- cosh 2 g H ( sinh q + V 



h(Af + 1)C sinh q-Aa 2 - A(Af - 2)a + Af 2 + 6Af - 4] 

8 cosh g 



- 2(a+jV-l) 



_/(Csinhg;a + A0 
where Vq is an irrelevant constant given by 



(sinhg-a + 3 1 . 2( sinh g + A/" + l 

2(a+Af-l 



/(Csinhg;a + A^) 2 _ 



4a 2 + 4(AT-4)a + AT 2 + 16 
Vo = 5 Co- 



(3.46) 
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Solvable sectors: 



e -C(sinhg)/2-Cgdq 

V M ~ = (ft (C smh g; a), . . . , ^(C sinh q; a)) ^ g)Af/2 _ 1/(c sinh ~ ~y (3.47a) 



Vat = <Xi(Csinhg;a+A/'),...,XAr(Csmhg;a + A/')) 

e C(sinhg)/2+Cgdg 

X (coshg)^/ 2 /(Csinhg;a + j\0' (3 ' 47b) 

where gdg = tan -1 (sinh q) is the Gudermann function. The solvability condition is not 
satisfied in this case and both of the Hamiltonians are only quasi-solvable. For a > 1, the 
potentials V ± (q) given in (13.451) are defined on the whole real line R. From the solvable 
sectors (JS3ZD, it is clear that neither VjJ )_ (M) nor V$ )+ (M) belongs to L 2 (R), so the JV-fold 
SUSY is dynamically broken in the constant-mass background. 

Now, the PDM potentials U ± {q) can be obtained with help of Eqs. (12.18a|) . f !3.45j) . and 
(I3.46p . and the solvable sectors Vjy of the corresponding PDM Hamiltonians H ± for an 
arbitrary mass function m(q) read from (j2.18cj) as 

Vat = (0i((smhu(q);a),...,<p A r((smhu(q);a) 



m(q) V4 e -C(sinh u(g))/2-C gd u(q) 

X (cosh u{q)yi 2 ~ l f{C, sinh u{q); a) ' (3 ' 48a) 
Y& = <Xi(Csinhu(g); a + Af), £aKC sinh «(<?); a + Af)) 
m(q) 1/4 e^ (sinh n (s))/ 2 +C § d u (g) 

X (cosh u{q))^/ 2 f{( smhu{q); a + Af)' (3 ' 48b) 

where u(q) is given by (I2.19p . Let us now consider two cases: 

(i) m(q) = sech 2 g, q £ (— oo, oo), for which the change of variable is u(q) = gdg. Then, 
the solvable sectors of t/ ± (g) are given by 

Va? = (^i(Csinhw(g);a),...,^(Csinhw(g);a;)) 
v^sech q e _1 > sinh (§ d i)/ 2 ~C gd(gd q) 



[cosh(gd q)}^ '/ 2 ~ l j '(( sinhw(g); a) ' 
V Af = (xi(Csinhw(g); a + Af),..., sinh u(q); a + Af)) 



(3.49a) 



\/sech g e^ sinh (§ d <2)/ 2 +C gd(gd g) 

X [cosh(gd q)]*/ 2 f{( sinh u(g); a+Af)' (3 ' 49b) 

In this case, the mass function as well as the potentials U ± (q) are well behaved on (— oo, oo), 
so we can consider the whole real line R as a domain of the potentials. From the solvable 
sectors (13 .49 p . it is clear that both Vjy-(R) belong to L 2 (R), which means unbroken Af- fold 
SUSY, i.e., the mass profile affects symmetry restoration. 

(ii) m(q) = 2e~ 2q /it. In this case, the solvable sectors Vjy reduce to 

Va> = <^i(Csinhu(g);a),. . . , ^(Csinhu(g); a)) 

exp[-g 2 /4 - Csinh(Erf q)/2 - Cgd(Erf q)} 

[cosh(Erfg)]^/ 2 - 1 /(Csinh«(g);«) ' 1 1 

V tf = <Xi(Csinhu(g); a + Af),..., XaKC sinh «(<?); a + Af)) 
exp[-g 2 /4 + Csinh(Erf q)/2 + (gd(Erf q)] 

[cosh(Erf q)]^/ 2 f(( sinh u(q); a + Af) ' 1 J 
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From the above solvable sectors (13.501) . it is clear that both V_^(M) belong to L 2 (R). That 
is, in this case we again have unbroken TV-fold SUSY. 

We note that there are other mass functions, e.g., m(q) = (j3 + q 2 ) 2 /(l + q 2 ) 2 , which have 
no effect on the dynamical breaking of A/"- fold SUSY, i.e., it is also possible to construct 
PDM systems which maintain the broken TV-fold SUSY. 



Example 3.6. A(z) = (z + C) 2 /2, [a 2 = 1/2, a x = ( = y/(a - l)(a + Af - 1) 
Potentials: 



V^-(q) 



l e 2 3 _ A^+l c , _ (jV-l)(jV F + 2a-2C-l)C c - q 



4 



+ 



( 2 [N 2 + 2TV(4a - 2( - 3) + 4a(2a - 2( - 3) + 4( + 5] 



-2q 



-2 



(a - C - l)e 9 2(a - l)e 2g 



/(e9-C;a) /(e*-C;a) 



(3.51) 



y(o)+(g) = ±e 29 + 



l„ 2q , A/--1 q , (A/ + l)(A/+2a-2C-l)C. 



-e y + 



+ 



( 2 [Af 2 + 2A/"(4a - 2C - 3) + 4a(2a - 2( - 3) + 4( + 5] 



-2g 



- 2 



(a + TV - C - l)e 9 2 (a + M - l)e 2q 



/(e«-C;a + AT) f^-C-a + N) 2 
where Vq is an irrelevant constant given by 



+ V , 



_ (TV + 2a) 2 + 2((N - 2a) + 2(7( - 8a + 8) 
~ 8 C °' 



Solvable sectors: 



V 



(0)- 



(<Pi{e q - C; a), • • • , ^(e 9 - ")) 



AT 



x exp 



/(e«-C;a) 
£ + 2C-2a-AT+l 
2 2 S 



TV-2 



■5 



(3.52) 



(3.53) 



V 



(o)+ _ (xi(e 9 -C;« + AT),...,XAr(e 9 -C;« + AA)) 



N 



x exp 



/(e«-C;a + A0 
e 9 2C-2a-A/"+l A _ 
Y 2 



Ce q - —q 



(3.54) 



This system is new and presented in this paper for the first time. The exponential-type 
Vjf(q) are naturally defined on the whole real line M since they have no singularity on it, so 
the Hilbert space is L 2 (M). Noting that 2( - 2a - TV + 1 < for a > 1, since 



4C 2 - (2a + TV - l) 2 = -4a - (TV - 1)(TV + 3) < -(TV + l) 2 < 0, 



we see that V_SJ )_ (M) C L 2 (R) and V$ )H 



(jL L 2 (R) for ( > 0. Hence, it manifests unbroken 



TV-fold SUSY. For C < 0, on the other hand, neither V$ )- (M) nor VjJ )+ (M) belongs to L 2 { 
so the TV-fold SUSY is broken in the constant-mass background. 
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In a PDM background, the solvable sectors Vj^ of the type X 2 PDM Hamiltonians H 
are deformed as [cf., Eq. f l2.18cj) 



± 



V 



AT 



{^(e^-Cia),...,^-^)) 
/(e«(*)-C;a) 

x m(g) 1//4 exp 



— — + — (e u(q) — u(q) 



(3.55) 



V 



<Xi(e u(9) - C; a + AO, • • • , XAKe u(9) - C; a + A/")) 



AT 



x m(g) 1//4 exp 



3"^) 2C-2a-A/"+l . A/" , / 

~2 2 Ce - Y uiq) 



(3.56) 



and the potentials U ± (q) can be obtained using Eqs. (|2.18ap . f)3.5ip . and (13.521) . In this 
case, the choice of mass functions are as follows: 

(i) m(q) = (1— q 2 ) , q € (—1, 1), for which the solvable sectors of the PDM Hamiltonians 
if ± are given by 



V 



x exp 



;i -g 2 ) 1/4 /(e M(9) -C;«) 



3 sin g 



+ 



2C-2a-A^+l 



sin q 



Af-2 



sin g 



(3.57) 



<Xi(e u(<?) - C; « + AO, • • • , xu^ u[q) - C; a + AT)) 



(1 -g 2 ) 1/4 /(e u(9) -C;« + A^) 



x exp 



2C-2a-AT+l _ sin -i a M . _ x 

— — 2 Ce y sm g 



(3.58) 



From the above solvable sectors, it is clear that both Vjy(— 1, 1) do not belong to L (— 1, 1), 
so it manifests broken TV-fold SUSY irrespective of the sign of (. Hence, comparing the 
normalizability conditions in both the constant and position-dependent mass cases, we con- 
clude that the mass function m(q) = (1 — q 2 )" 1 affects dynamical breaking of TV-fold SUSY 
for C > 0. 

(ii) m(q) = 2e~ 2q /tt, q E (— oo, oo), for which the TV-fold SUSY remains unbroken, which 
is evident from the corresponding solvable sectors given by 



V 



(ft(e«fo>-C; «)>■■■> ^(e n(9) -C;a)> 



AT 



/(e«<9)-C;a) 



exp 



' 4 



+ 



(3.59) 



V 



A r 



<Xi(e u(9) - C; a + AT), . . . , XAr(e"^ - C; a + AO) 



/(e«(ff) -C;« + A0 



x exp 



g 2 e Erf9 2C - 2a - 
"T + "1 2~ 



^Ce^-^Erf, 



(3.60) 
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From the normalizability conditions in the constant and position-dependent mass cases, we 
see that the mass function m(q) = 2e~ 2q /it affects the dynamical breaking of A^-fold SUSY 
for C < 0. 

IV. SUMMARY AND PERSPECTIVES 

In this paper, we have investigated effect of position-dependent mass background on 
dynamical breaking of type B and type X 2 TV-fold SUSY. We have selected three different 
models in the constant mass background for each type, and then examined whether some of 
the physically relevant effective mass profiles can affect the pattern of A/"-fold SUSY breaking 
in each model. We summarize the results in Table [H We can easily see from Table U that, 
except for the rational potentials, some of the PDM profiles can actually affect and change 
the patterns of dynamical /V-fold SUSY breaking in all the trigonometric, hyperbolic, and 
exponential potentials. Although we have selected the specific types of A/"-fold SUSY to 
develop physical applicability of the new mathematical concept of exceptional polynomial 
subspaces, we can of course make a similar analysis on other types of A^-fold SUSY such as 
type A and type C to find out positive effect of PDM on SUSY breaking in some models. 

Hence, it would be possible to observe experimentally transition between a broken and an 
unbroken phases if an effective mass can be controlled experimentally such that the constant 
mass limit can be also realized in an experimental setting. The physical meanings of a 
position-dependent mass depend on each physical system under consideration, for instance, 
the curvature of the local band structure of an alloy in the momentum space for electrons 
in a crystal with graded composition the particle densities of 3He and 4He in pure and 
mixed helium clusters with doping atoms or molecules 0], the effective electron mass for 
electrons confined in a quantum dot [IJ and for dipole excitations of sodium clusters jEf, and 
so on. Thus, if we can prepare such an atomic, molecular, or condensed matter system which 
is described by a certain PDM quantum model subjected to an A/"-fold SUSY potential with 
mass profiles, e.g., m(q) = e~ u q or (1 — v 2 q 2 )~ l where v is an experimentally adjustable 
parameter such that v — > is realizable, then the spectral change of the system could be 
observed at v — due to the phase transition. The essence and novelty of our idea rely 
on the observation that the physically controllable PDM can cause the phase transition 
by changing the normalizability of the solvable sector although the latter is superficially a 
simple mathematical aspect. Hence, it is quite important to note that the normalizability 
of wave functions can play much more roles than the quantization of energy spectrum which 
is referred to by any standard textbook on quantum mechanics. 

We note that this experimental observability might have impact not only on some atomic, 
molecular, and condensed matter problems from which PDM quantum theory originated, but 
also on high-energy physics. Until now many high-energy physicists have believed that SUSY 
is realized at the GUT or Planck scale as a resolution of the naturalness and the hierarchy 
problem but is broken at least at the electroweak scale. Unfortunately, however, theoretical 
analysis on dynamical SUSY breaking in field theoretical models are extremely difficult on 
the one hand, and it is virtually impossible to make a GUT scale experiment on the other 
hand. The aforementioned experimental observability suggests that we might extract some 
clues to understand dynamical SUSY breaking in high-energy physics from realistic eV scale 
experiments in atomic, molecular, and condensed matter physics. It is because the Witten's 
work [lH has indicated that the mechanism of dynamical SUSY breaking in quantum field 
theory and quantum mechanics is essentially the same. We also note that the careful non- 
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TABLE I: The effects of PDM profiles on dynamical breaking of TV-fold SUSY in various type B 
and type X2 models. 



Types of potentials 


Dynamical breaking of jV-fold SUSY 






Constant mass 


PDM 




rational 


unbroken 


no effect 


Type B 


trigonometric 


broken 


unbroken for m(q) oc e~ 2qi 






partially broken for 


broken for m(q) = (1 — q 2 )^ 1 




exponential 


-M/2<b 1 < (M + l-2k)/2 
and fully broken otherwise 


unbroken for m(q) oc e~ 2q2 




rational 


unbroken 


no effect 


Type X 2 


hyperbolic 


broken 


unbroken for m(q) oc sech 2 q, e~ 2qA 




exponential 


unbroken for £ > 
broken for £ < 


broken for m(q) = (1 — q 2 ) , V£ 
unbroken for m{q) oc e~ 2q2 , V£ 



perturbative analyses in Refs. 59|, |60| have shown that the mechanism of dynamical breaking 
of ordinary and TV-fold SUSY is also the same. Hence, dynamical aspects of SUSY quantum 
field theoretical models would be mimicked in TV-fold SUSY quantum mechanical toy models, 
regardless of whether or not TV-fold SUSY can be realized in higher dimensions. Therefore, 
we believe that further studies in this direction are worth pursuing both theoretically and 
experimentally. From a theoretical point of view it is a challenging issue to investigate both 
a perturbation theory and the non-renormalization theorem in PDM quantum systems. 
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